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Abstract 

We give the correct analytic expression of a finite integral appearing in the four-loop computation 
of the renormalization-group functions for the two-dimensional nonlinear cr-model on the square lat- 
tice with standard action, explaining the origin of a numerical discrepancy. We revise the numerical 
expressions of Caracciolo and Pelissetto for the perturbative corrections of the susceptibility and of 
the correlation length. For the values used in Monte Carlo simulations, N = 3, 4, 8, the second 
perturbative correction coefficient of the correlation length varies by 3%, 4%, 3% respectively. Other 
quantities vary similarly. 
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The two-dimensional nonlinear u-model has been extensively studied since it is supposed to share with four- 
dimensional QCD the property of being asymptotically free. A lot of work has been devoted to the checking of 
renormalization-group predictions. In order to make accurate comparisons it is very important to have high-order 
perturbative predictions. In M the renormalization-group (3- and 7-functions were computed for the standard action 
with nearest-neighbour interactions on a square lattice. The result was expressed in terms of a small number of 
basic lattice integrals that were computed numerically. In Ref. M, using the coordinate space method 0, the lattice 
integrals that appear in the results of |jj were independently evaluated with higher numerical precision, and two of 
them were found grossly incorrect. As we shall discuss below, the large discrepancy found for one of them was due to 
the fact that its definition was not correctly given in Ref. jjj . 

Two of the present authors (B. A. and M. P.) have performed a complete independent check of the analytic part 
of the lattice calculation, by computing the two-point function in the presence of an external constant magnetic field 
to three loops. Their result is in agreement with that presented in Ref. fjj. However their computation allowed the 
discovery of an incorrect definition in Ref. jjj : the constant W2 which was used in the computation of Ref. M does 
not have the definition given in App. A.l. Let us indicate with W 2 the definition appearing in App. A.l of Ref. jjj: 
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where p 2 = ~^2 4sin 2 p^/2, 1(h) is the integral of the propagator 
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and R w 0.0148430 is a numerical constant introduced in m. In the calculation of Ref. fil, W2 was instead defined as 
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This difference expiains most of the discrepancy found in Ref. H for the numerical estimate of W%. Using 0] 

W 2 = 0.0006923019(1), (0.4) 

wc find 



W 2 = 0.002122552, 



(0.5) 



which should be compared with W 2 = 0.00221 reported in Ref. M. The remaining small discrepancy was due to 
numerical problems that were not understood at the time. We have completely revised the programs and now we find 
W% = 0.002122(1) in agreement with (0.5). It follows that, while the analytic expressions of Ref. O are correct, the 
numerical results for the perturbative corrections presented in Refs. |l|,[2| must be revised. Below we give the numerical 



expressions for the perturbat ive c onstants a 2 , b 2 , C3, and b 2 
of Ref. and the estimate (0.5): 
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(see Ref. M for definitions), using the numerical estimates 
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(0.0444 + 0.0216iV + 0.0045iV 2 - 0.01297V 3 ), 
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For N — 3 we have a 2 = —0.1969, b 2 — —0.2853, C3 = 0.1346: with respect to the previous estimates of Ref. [[JJ , a 2 , 
b 2 , and C3 vary by 3%, 5%, and 1% respectively. The difference is very small and does not change the conclusions of 
the papers that compared Monte Carlo results with the perturbative predictions p-fl . 
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